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ABSTRACT 

There is increasing evidence for the presence of blue super star clusters in the 
central regions of ultraluminous infrared galaxies like Arp 220. Ultraluminous 
galaxies are thought to be triggered by galaxy mergers, and it has often been 
argued that these super star clusters may form during violent collisions between 
gas clouds in the final phase of the mergers. We now investigate another set of 
models which differ from previous ones in that the formation of the super star 
clusters is linked directly to the very intense starburst occurring at the very 
center of the galaxy. Firstly we show that a scenario in which the super star 
clusters form in material compressed by shock waves originating from the central 
starburst is implausible because the objects so produced are much smaller than 
the observed star clusters in Arp 220. We then investigate a scenario (based on 
the Shlosman-Noguchi model) in which the infalling dense gas disk is unstable 
gravitationally and collapses to form massive gaseous clumps. Since these 
clumps are exposed to the external high pressure driven by the superwind (a 
blast wave driven by a collective effect of a large number of supernovae in the 
very core of the galaxy), they can collapse and then massive star formation 
may be induced in them. The objects produced in this kind of collapse have 
properties consistent with those of the observed super star clusters in the center 
of Arp 220. 

Subject headings: galaxies: individual (Arp 220) - galaxies: starburst - stars: 
formation 
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1. INTRODUCTION 

Galaxy mergers provide a mechanism for the metamorphosis of galaxies and sometimes 
trigger intense star formation events (starbursts) in the central regions of the merger 
remnants (Toomre 1977; Schweizer 1982). When such intense starbursts happen, 
ultraluminous infrared galaxies (ULIGs) result because of the enormous energy released 
from the central regions of the galaxies (Soifer et al. 1984; Wright et al. 1984; Sanders et 
al. 1988; Sanders & Mirabel 1996). One of the remarkable properties found in such galaxy 
mergers is that a number of blue super star clusters (SSCs) are resided in their central 
regions (Lutz 1991; Ashman & Zepf 1992; Holtzman et al. 1992; Shaya et al. 1994; see 
for a review Zepf & Ashman 1993). One plausible idea is that such SSCs may be formed 
by violent collisions among gas clouds in the final phase of the galaxy mergers (Ashman & 
Zepf 1992; Zepf & Ashman 1993). 

The SSCs in the mergers are then distributed over the circumnuclear regions which 
are spatially close to the currently active starburst core (cf. Lutz 1991; Shaya et al. 1994). 
Such SSCs, while common in the cores of ULIGs, are rarely seen in galaxies where there 
is no evidence for an intense burst of star formation at the very center. Therefore, it is 
intriguing to investigate physical relationships between the circumnuclear SSCs and the 
nuclear starburst and address the question: is the formation of SSCs attributed to some 
effects of the nuclear starburst? This is the subject of the present paper. The link between 
the concentrated starburst in the core and the formation of the SSCs is a fundamental 
difference between the models we discuss and those that rely on violent collisions mentioned 
in the previous paragraph. 

The nearest and best-studied ULIG is Arp 220 (z = 0.02) and we use it as a basis for 
comparison in the present work. In this galaxy, two nuclei exist with a projected separation 
of 350 pc (e.g., Condon et al. 1991). Recent optical and near-infrared imaging made with 
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the Hubble Space Telescope has shown that about a dozen of star forming clumps are 
distributed in the circumnuclear region with a radius of several hundred pc in Arp 220 
(Shaya et al. 1994; Scoville et al. 1998). In particular, the western nucleus of Arp 220 is 
surrounded by several clumps which are distributed along the north-south direction. In 
Figure 1 we show a schematic illustration of the spatial distribution of SSCs in Arp 220. 
These clusters are distributed within a radius of 700 pc from this western nucleus. The 
masses of SSCs are estimated to be ~ 10 6 - 10 8 M Q although these numbers depend on the 
assumed stellar populations (Shaya et al. 1994). 

These SSCs are more massive and larger (~ 10 pc) than typical globular clusters in 
the Galaxy. Nevertheless, it is possible that the formation of SSCs provides important 
hints on the formation of globular clusters. Recent theoretical studies have suggested that 
globular clusters are envisaged as forming in giant molecular clouds, in the same way as 
we see star-cluster formation in our Galactic disk (Harris & Pudritz 1994; McLaughlin & 
Pudritz 1996; Elmegreen & Efremov 1997). Early in the formation of the Milky Way, which 
is when the globular cluster formed, there may also have been considerable concentrated 
star formation, as the bulge and the central part of the disk formed. This hypothetical 
comparison provides a secondary motivation for investigating the possible role of these 
concentrated starbursts in galaxy cores in forming massive star clusters. 

2. STARBURST-DRIVEN STARBURSTS IN ARP 220 

Arp 220 shows a double-bubble structure of extent 13 kpc (Heckman et al. 1987, 1996 
see Figure 1). Such a morphology is evidence for a superwind i.e. a the blast wave driven by 
a large number of supernovae (Tomisaka & Ikeuchi 1988; Heckman, Armus, & Miley 1990; 
Suchkov et al. 1994). A typical galactic-scale nuclear starburst forms numerous (> 10 4 ) 
massive stars during a short period (~ 10 7 years; Weedman et al. 1981, Balzano 1983). 
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Therefore, a burst of supernovae occurs ~ 10 7 years after the onset of the starburst. Since 
these numerous supernovae release a huge amount of kinetic energy into the circumnuclear 
gas, the circumnuclear gas is thermalized and then ejected from the starburst region as a 
superwind. This superwind interacts with ambient gaseous matter, leading to the formation 
of shocked gaseous shell in the circumnuclear region. If there is a dense gaseous disk 
surrounding the nuclear starburst regions, this disk can collimate the superwind along the 
direction perpendicular to the disk (Tomisaka & Ikeuchi 1988); such a disk is observed in 
Arp 220 by Scoville et al. (1997). Gas in the inner parts of this dense collimating disk is 
will be affected by both the dynamical (e.g., shocks) and thermal (e.g., compression by the 
high-pressure hot gas) effects of the superwind and it is these processes that we consider in 
the rest of this paper. The physical processes we describe here are quite distinct from what 
might be happening in the outer regions of the galaxy where the superwind interacts with 
the ambient hot diffuse coronal gas. At this stage we note that in Figure 1, the SSCs are 
distributed along the direction almost perpendicular to the double-bubble structures; this 
is what we expect if the SSCs form in the dense collimating gas disk. 

The inner parts of the dense gas disk should experience a significant dynamical effect 
from the superwind. While direct evidence of shocked gas in Arp 220 is not observed, 
we note that shocked regions are indeed observed in the nearer but less luminous nuclear 
starburst galaxy M82 (Lester et al. 1990). Although the actual gaseous matter in the nuclear 
region of Arp 220 is almost certainly clumpy, it is instructive to derive order-of-magnitude 
estimates of the dynamical effects of the superwind by considering the effects of a superwind 
propagating through a homogenous medium. The supernovae responsible for shocking the 
gas occur continuously over a timescale longer than or comparable with the dynamical 
timescale of the initial gas cloud, and the evolution of the shocked material can be described 
by a superbubble model (McCray & Snow 1979; Koo & McKee 1992a, 1992b; Heckman et 
al. 1996; Shull 1995 and references therein). The radius and velocity of the shocked shells 



- 6- 



formed in the dense gas disk at time t are then 



r she n - 440L m / ech)43 7i H;4 / t/ pc, (1) 

and 



v shcll ~ 26L 1 J l ^ chA3 n H J 5 t 7 2/5 km s \ (2) 

where L m ech,43 is the mechanical luminosity released collectively from the supernovae in 
units of 10 43 erg s" 1 , nn,4 is the average hydrogen number density of the ISM in the dense 
gas disk, assumed uniform, in units of 10 4 cm -3 , and t 7 is the elapsed time since the onset 
of the superwind in units of 10 7 years. For Arp 220, recent observations suggest that 
£ m cch,43 = 1 (Heckman et al. 1996) and n H ,4 = 1 (Scoville et al. 1991, 1997). Although we 
do not know physical conditions of nuclear gas disk in the recent past (i.e., ~ 10 7 years 
ago), they may not be dramatically different from those observed currently. We therefore 
assume riu = 10 4 cm -3 . In any case, the dependencies on riu of both r s h e ii and i> s h e ii are 
very weak (i.e., oc n^ 5 ). Hence if our assumed values of wh is incorrect by a factor of two 
or three, the effects on the shell properties are very weak. In order to explain the spatial 
extension of the double-bubble structures, Heckman et al. (1996) estimated an age of the 
superwind, t 7 ~ 3 given a low-density ambient gas (% ~ 10~ 2 cm -3 ). Using this age, we 
obtain r shell ~ 800 pc for the dense gas disk. This is approximately the extent of the SSCs. 

Next we investigate the formation mechanism of the SSCs in the dense gas disk around 
the western nucleus. Gaseous fragments can grow gravitationally in the shocked shell formed 
by the dynamical effect of the superwind (cf. Whitworth et al. 1994). Fragments that form 
within the shells experience a net inward acceleration due to self-gravity and a net outward 
acceleration due to the internal pressure. Whitworth et al. (1994) investigate the balance 
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between these two accelerations and show that the time-scale for the growth of the fastest- 
growing fragments is tf as tcst ~ (2c s )/(G£) where c s is the sound speed in the shell, and £ is 
its surface mass density. Non-linear fragmentation in the shell then happens first at a time 
t = tf as test- Noting that the surface mass density £ = CnH^H^sheii, where C is a constant 
determined by the geometry (C = 1/3 for a sphere) and m H is the hydrogen atom mass, we 
then find from equations (1) and (2) that fragmentation within the shell first happens at a 
time t { ~ 5.7 x l^C^ 2 n^[ 2 M'^ 2 yr at a radius r f ~ SOL^ ^Co ^n^Mifl 10 pc. 
Here Mf : io is the Mach number when these fragments first appear in units of 10, equal to 
^sheii/cs (here it is assumed that M{ 3> 1; Whitworth et al. 1994). Estimating c s is difficult 
because the turbulent pressure in the shell is much greater than the thermal pressure. From 
equation (2), the shell is moving outward at a velocity V{ when the fragments first appear, 

1/5 1/5 1 

where Vf ~ 81 43 A4 f 10 km s . Thus V{ is almost independent of c s , so that our lack 
of knowledge of the sound speed is unimportant in determining the velocity of the shell 
when the fragments form. Following Whitworth et al. (1994), we can estimate the mass 
and size of the fragments: M f ~ cl^ 2 (G 3 n-nm} 1 Vi)~ 1 ^ 2 ~ 7.9 x l^cj^n^^vj^l 2 M© and 
/f ~ c^ 2 (Gnnm-ftVf)' 1 / 2 ~ ^■ / ^c^ 2 n^ i 1 { 2 v^g{ 2 pc where c s and t>f are expressed in units of 10 
km s" 1 and 81 km s _1 , respectively. We note that the local sound speed in the molecular 
gas clouds is c' ocal ~ 0.5T^ 2 km s~ x where T 50 is the temperature of molecular gas in units 
of 50 K. However, when we discuss the fragmentation in the shell, it is more reasonable 
to regard the random velocity of clouds as an effective sound speed - this bulk velocity 
is much larger than the local thermal velocity (e.g., 90 km s _1 in Arp 220; Scoville et al. 
1998). The conclusion from this simple calculation is therefore: fragments of gas that grow 
gravitationally in the shocked shell are too small (~ 10 4 M© and ~ 5 pc) to be the SSCs 
observed in Arp 220. Adopting a smaller sound speed (which is the main uncertainty in the 
above calculation) makes this statement even stronger. 

Therefore we have to look for another formation mechanism of the progenitor of SSCs. 
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We consider the instability of the dense gas disk in the central region of Arp 220 due to its 
self gravity. We draw on the results of Shlosman & Noguchi, who investigated the general 
problem of the gravitational instability of nuclear gas disks in detail. They estimated the 
mass of a superclouds forming in such a disk as 

M scl = vr^ 2 S g (3) 

where S g is the surface mass density of the gas disk and A cr i t is a typical scale length of the 
supercloud given by 

A„„ = ^ (4) 

K 

where k is the epicyclic frequency at the region where the gravitational instability occurs. 
Equations (3) and (4) can be expressed: 



M scl ~ 5.7 x lO 9 !]^ 4 M , (5) 

and 



A crit ~ 8.5 x 10 2 E g;4 «;- 2 pc, (6) 

where S gj4 is in units of 1O 4 M pc~ 2 and k 3 is in units of 10 3 km s" 1 kpc -1 . Using the 
observed values of S g and k (Scoville et al. 1998), we obtain masses of the superclouds up 
to ~ 1O 9 M (see Table 1). These are the masses of the gaseous clouds which will only equal 
the mass of forming stars, if the star formation efficiency is 1. Adopting a more reasonable 
efficiency of 0.1 from observations of the Galaxy, we derive a final stellar mass of ~ 1O 8 M , 
which is consistent with observation (Shaya et al. 1994). As shown in Table 1, the masses 
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of the gaseous clouds are less massive with increasing with radius. On going from 0.1 kpc 
to 1 kpc the clusters are approximately an order of magnitude smaller assuming the same 
star formation efficiency everywhere. Despite the (large) uncertainty due to extinction, this 
trend is close to what is seen by Shaya et al. (1994) - see their Tables 1 and 3. 

Finally, we consider how the massive and large proto-SSC clouds that came from our 
stability analysis evolve into the more compact and less massive SSCs seen in the HST 
images. Once the superclouds are formed, they will collapse because of the effect of external 
pressure (Kimura & Tosa 1993). Given the hydrogen number density and the kinetic 
temperature of the superclouds, % ~ 10 4 cm -3 and T^n ~ 10 K, the pressure supported by 
the random motion of gas is P mt ~2x 10~ n dyn cm -2 . Next we estimate the pressure of 
inter-supercloud gas. Since the hot gas traced by soft X-ray emission is indeed associated 
with the nuclear region of Arp 220, it seems reasonable to assume T kin ~ 10 7 K (Heckman 
et al. 1996). However, there is no direct information on the inter-supercloud gas density 
^h,ic- Although there is considerable amounts of gas (~ 1O 1O M ) in the nuclear region, 
most of this gas occurs in dense molecular clouds (Scoville et al. 1997). The only statement 
that we can make with any confidence about the density of intercloud gas is that it must 
be intermediate between the cloud density (~ 10 4 cm~ 3 ) and the density of the ambient 
coronal gas (~ 10~ 2 cm -3 ; Heckman et al. 1996). A reasonable value is probably around 1 
cm -3 , which is the typical gas density of the cold ISM in our galaxy (Spitzer 1978). The 
external pressure provided by the hot gas is then P cxt ~2x 10" 9 (nH,ic/lcuT 3 ) dyn cm -2 
for T kin ~ 10 7 K. Therefore, the external pressure caused by the hot gas is higher than the 
internal one for any value of the intercloud density higher than that of the coronal gas, 
and is much higher (by a factor of 100) if n^ic ~ 1cm -3 . Since the cooling timescale of 
neutral gas in the fragments, t coo i ~ 100(n H /10 4 cm" 3 )" 1 yr (Spitzer 1978), is quite short, 
stars could form just after the fragmentation. The stars may form through the process of 
gravitational instability occurring inside the fragment. A typical length of the gravitational 
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fluctuation (the Jeans' length: Jeans 1929) is estimated to be Aj ~ (^^""V^Pgas) 1 ^ 2 ~ 

1 pc where p gas = numu. Therefore, since a typical mass of the fluctuation amounts 
torn ~ -^jPgas ~ 1OOM , it is likely that the fluctuation will evolve into a massive star 
with a mass of ~ 1OM . The lifetime of these massive stars is ~ 10 7 years. The single 
fragment (comprising these stars) experiences dynamical relaxation on a time scale of 

_ 3/2 1/2 

2 x 10 AJ it 100 M scl g years where A critj ioo is the diameter of the supercloud in units of 100 
pc and M sc i j8 is the mass of the supercloud in units of 1O 8 M . These predicted values of the 
cloud size and mass are close to the observed ones (several 10 pc and 10 6 - 1O 8 M ; Shaya 
et al. 1994) in Arp 220. 

Therefore adopting a Shlosman-Noguchi model to investigate the gravitational 
instability of the circumnuclear disk in Arp 220 allows us to explain the formation of SSCs 
found there. The link between the formation of the clusters and the central starburst that 
was discussed in Section 1 is clear - the supernovae occurring in the center lead to the 
formation of the superbubble, which in turn provides an external pressure large enough to 
initiate collapse in the gas disk. 
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Figure captions 



Fig. 1: The blue stellar clusters found in Shaya et al. (1994) are overlaid on the 4.83 
GHz radio continuum map obtained by Baan and Haschick (1995). The blue clusters are 
mostly distributed along the NS direction through the western radio nucleus. A pair of the 
large-scale emission-line bubbles, traced by Ha + [Nil] A6548, 6584 emission, is shown in the 
lower panel taken from Heckman et al. (1996). The bubbles are roughly elongated along 
the EW direction which is perpendicular to that of the blue clusters. 



Table 1: Physical properties of the superclouds formed in the dense gas disk 



R 


S g 


K 


Acrit 


M scl 


(kpc) 




(km s _1 kpc -1 ) 


(P C ) 


(M Q ) 


0.1 


3.5 x 10 4 


3.6 x 10 3 


230 


1.45 x 10 9 


0.3 


0.8 x 10 4 


1.5 x 10 3 


302 


5.76 x 10 8 


0.5 


0.4 x 10 4 


1.0 x 10 3 


340 


3.65 x 10 8 


1.0 


0.1 x 10 4 


0.5 x 10 3 


340 


9.12 x 10 7 




90 arc sec 



